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Branched glycerol dialkyi glycerol tetraethers (bGDGTs) are membrane-spanning lipids 
that likely stabilize membranes of some bacteria. Although bGDGTs have been reported 
previously in certain geothermal environments, it has been suggested that they may 
derive from surrounding soils since bGDGTs are known to be produced by soil bacteria. 
To test the hypothesis that bGDGTs can be produced by thermophiles in geothermal 
environments, we examined the distribution and abundance of bGDGTs, along with 
extensive geochemical data, in 40 sediment and mat samples collected from geothermal 
systems in the U.S. Great Basin (temperature: 31-95°C; pH: 6.8-10.7). bGDGTs were 
found in 38 out of 40 samples at concentrations up to 824ng/g sample dry mass and 
comprised up to 99.5% of total GDGTs (branched plus isoprenoidal). The wide distribution 
of bGDGTs in hot springs, strong correlation between core and polar lipid abundances, 
distinctness of bGDGT profiles compared to nearby soils, and higher concentration of 
bGDGTs in hot springs compared to nearby soils provided evidence of in situ production, 
particularly for the minimally methylated bGDGTs I, lb, and Ic. Polar bGDGTs were found 
almost exclusively in samples <70°C and the absolute abundance of polar bGDGTs 
correlated negatively with properties of chemically reduced, high temperature spring 
sources (temperature, H2S/HS^) and positively with properties of oxygenated, low 
temperature sites (O2, NO3 ). Two-way cluster analysis and nonmetric multidimensional 
scaling based on relative abundance of polar bGDGTs supported these relationships and 
showed a negative relationship between the degree of methylation and temperature, 
suggesting a higher abundance for minimally methylated bGDGTs at high temperature. 
This study presents evidence of the widespread production of bGDGTs in mats and 
sediments of natural geothermal springs in the U.S. Great Basin, especially in oxygenated, 
low-temperature sites (<70°C). 
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INTRODUCTION 

Glycerol dialkyi glycerol tetraethers (GDGTs) are core membrane- 
spanning lipids that resist delamination and impart membrane 
stability at high temperature and low pH. GDGTs are synthe- 
sized by a variety of extremophiles and some non-extremophiles 
and occur widely in nature (Schouten et al, 2000, 2013 and 
references therein). The isoprenoid GDGTs (iGDGTs) contain 
alkyl groups constructed by polymerization of isoprene sub- 
units and are synthesized by physiologically and phylogenetically 
diverse archaea [reviewed in Schouten et al. (2013)]. A sepa- 
rate group of GDGTs, the branched GDGTs (bGDGTs), have 
been identified in lipid extracts prepared from a wide variety 
of environments, including peat bogs, soils, estuaries, and lake 
and river water and sediments (Schouten et al., 2000; Sinninghe 



Damste et al., 2000, 2009; Weijers et al, 2006a, 2009; Zink 
et al, 2010; Zhang et al., 2012; Liu et al, 2013; Yang et al, 
2011, 2013; Wang et al, 2013). bGDGT types vary according to 
the degree of methylation of the alkyl chains, with 4-6 methyl 
groups per GDGT, and the number of cyclopentyl moieties, with 
0 to 3 rings per GDGT (Figure SI) (Sinninghe Damste et al., 
2000). 

Although bGDGTs contain the membrane-spanning feature 
and ether linkages similar to archaeal iGDGTs, they are not 
isoprenoidal and differ in the stereochemical configuration of 
the second carbon position of the glycerol backbone; therefore, 
bGDGTs were proposed to have a bacterial rather than archaeal 
origin (Weijers et al., 2006b). Following a number of investi- 
gations focusing on microorganisms inhabiting water-saturated. 
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anaerobic peat bog environments (e.g., Weijers et al., 2009), a 
single bGDGT, bGDGT I (Figure SI), was discovered in pure cul- 
tures of Acidobacteria, specifically Edaphobacter aggregans Wbg- 
1^ and Acidobacteriaceae strain A2 - 4c (Sinninghe Damste et al., 
2011). The origin of the other bGDGTs found in peat bogs and 
other natural environments (Figure SI) remains unknown. Yet- 
uncultivated Acidobacteria are candidates for these other bGDGTs 
since all tested Acidobacteria produce abundant 13,16-dimethyl 
octacosanedioic acid (wo-diabolic acid), a likely precursor in 
bGDGT synthesis (Sinninghe Damste et al., 2011), and since var- 
ious bGDGTs and yet-uncultivated Acidobacteria co-exist in peat 
soils (Weijers et al, 2009). 

Although bGDGTs are known predominantly from soils, a few 
studies have documented bGDGTs in terrestrial geothermal envi- 
ronments. Substantial amounts of bGDGTs, accounting for up 
to 64% of total GDGTs, were recovered in all nine geothermal 
spring samples collected by Schouten et al. (2007) in Yellowstone 
National Park. However, the authors suggested that the majority 
of bGDGTs were from soil runoff (Schouten et al., 2007). A sub- 
sequent study of bGDGT lipids in soil transects near two springs 
in Surprise Valley, California, revealed the presence of bGDGTs 
in geothermally heated soils that were 12-41°C at the time of 
collection (Peterse et al., 2009). For both springs, bGDGT con- 
centrations in heated soils near the springs were greatly elevated 
over cooler, more distant soUs. bGDGTs were also recovered from 
two alkaline hot springs in Tibet, ranging from 52.0-83.6° C (He 
et al., 2012). The high concentrations of bGDGTs in some sam- 
ples and a unique lipid profile compared with nearby soils were 
consistent with bGDGT production in the springs. Most recently, 
a study of GDGTs extracted from natural sediments and incu- 
bation experiments in Great Boiling Spring (62-82° C), Nevada, 
provided multiple lines of evidence demonstrating the in situ pro- 
duction of bGDGTs in this hot spring (Zhang et al., 2013). In 
addition, cellulosic material that was incubated in situ in both 
sediments and the water column at 77 and 85° C also contained 
significant concentrations of bGDGTs. The bGDGT data were 
compared with 16S rRNA gene pyrotag datasets obtained from 
the same samples (Cole et al., 2013; Peacock et al., 2013), revealing 
that Acidobacteria were rare in all samples (<0.01% of pyrotags), 
and suggesting that the bGDGTs were produced by thermophilic 
bacteria, possibly yet-uncultivated Bacteroidetes, candidate phy- 
lum EM3 (Rinke et al., 2013), or candidate phylum "Atribacteria" 
(formerly candidate phylum OP9; Dodsworth et al., 2013; Rinke 
etal, 2013). 

In this study, we sought to augment the finding by Zhang 
et al. (2013) and to further test the hypothesis that bGDGTs 
are produced in natural geothermal environments and to exam- 
ine relationships between bGDGT composition and abundance 
and physicochemical setting. To do this, we examined the dis- 
tribution and abundance of bGDGTs in natural sediment and 
mat samples collected from a variety of geothermal springs in 
the northwest Great Basin (U.S.A.), along with a large physic- 
ochemical dataset. bGDGTs were widespread in the springs but 
were more abundant in cooler, more oxidized springs with 
well-developed microbial mats, defining the optimal habitat for 
bGDGT-producing thermophiles in habitats <70°C with abun- 
dant biomass. 



MATERIALS AND METHODS 

SAMPLE SITES AND PHYSICOCHEMICAL MEASUREMENTS 

Sediment and mat samples (top 1 cm of sediment/mat and 
water interface) were collected from eight hot springs located 
in northwest Nevada and northeast California (Figure 1). The 
same samples were used for simultaneous extraction of bGDGTs 
and iGDGTs; the results for iGDGTs are presented separately 
(Paraiso et al., in review). The sediment/mat water interface (top 
~1 cm) was collected. All samples were frozen on dry ice in the 
field and were transported and stored frozen (— 80°C) until they 
were thawed for lipid extraction and analysis. Field geochemical 
measurements and samples for laboratory geochemical measure- 
ments were taken prior to sediment and mat sampling (Costa 
et al, 2009; Vick et al., 2010). Briefly, temperature, pH, and con- 
ductivity were determined using a pH probe with temperature 
correction (LaMotte 5 Series, Chestertown, MD or YSI Model 
30, Yellow Springs, OH and WTW Model pH330i, Weilheim, 
Germany). Redox-sensitive analytes (O2 and sulfide) were mea- 
sured in the field using commercial kits (Hach, USA) with 
modifications for high temperature according to Miller- Coleman 
et al. (2012). Anions and cations were analyzed in the lab by ion 
chromatography (Dionex DX-500 chromatograph, AS14A col- 
umn, with 10 |xM Na2C03/NaHC03 as an eluent, Dionex, USA) 
and direct current plasma emission spectrometry (DCP-OES, 
Beckman, USA). Dissolved inorganic nitrogen species (NO3", 
NO2", NH3/NH+) were measured by automated colorimetry 
(Lachat, USA) as described by Dodsworth et al. (201 la,b). 

LIPID EXTRACTION AND LIQUID CHROMATOGRAPHY-MASS 
SPECTROMETRY (LC-MS) 

bGDGTs and iGDGTs were extracted from 5 gram lyophilized 
samples using a modified Bligh-Dyer extraction method as 
detailed by Lengger et al. (2012). The resulting lipids were sep- 
arated by silica-gel column chromatography using n-hexane:ethyl 




FIGURE 1 I The study area (inset, orange) was located with the U.S. 
Great Basin. Hot spring sampling sites (orange diamonds) were located 
within California and Nevada. 
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acetate (1:1, v:v) to collect the nonpolar fraction (Fl) and MeOH 
to collect the polar (F2) fraction. The polar fraction was divided. 
One aliquot, fraction F2A was not processed further. The other 
fraction, F2B was hydrolyzed with MeOH:HCl (95:5; v:v) at 70°C 
for 3h, extracted with DIH2O and DCM. All three fractions 
were dried under N2, dissolved in n-hexane:isopropyl alcohol 
(99:1; v:v), filtered through a 0.45 |xm polytetrafluoroethylene fil- 
ter, and dried under N2. Dried lipids were dissolved in 600 |xL 
of n-hexane:isopropyl alcohol (99:1; v:v) for analysis. Fl and 
F2A were run directly on LC-MS and F2B was run on LC-MS 
after hydrolysis. F2A was to check whether any core GDGTs may 
be present in the polar fraction. If there are any core-GDGTs 
detected in F2A, they will be subtracted from F2B and added 
toFl. 

Each fraction was spiked with a GDGT C46 internal stan- 
dard and analyzed on an Agilent 1200 liquid chromatography 
equipped with an automatic injector coupled to QQQ 6460 
MS and Mass Hunter LC-MS Manager software. Separation of 
GDGT peaks was achieved using a Prevail Cyano column (2.1 x 
150 mm, 3 (Jtm; AUtech, Deerfield, IL, USA) as described in Zhang 
et al. (2012) and GDGTs were quantified using Agilent 6460 
triple-quadrupole MS with an atmospheric pressure chemical 
ionization (APCI) source as (Zhang et al., 2012). 

STATISTICAL ANALYSES 

Two-way cluster analyses and non-metric multidimensional scal- 
ing (NMS) were completed in PC-ORD (MjM Software Design). 
Input data for these analyses included the relative abundance 
of the polar fractions of individual bGDGTs. Within PC-ORD, 
lipid data were normalized according to the maximum value of 
each lipid and were used to calculate a S0rensen (Bray-Curtis) 
distance matrix, which was applied in NMS and two-way clus- 
tering analyses. For two-way cluster analyses, the flexible beta 
method (P = —0.25) was used to create an agglomerative hier- 
archical clustering tree, with methods modified from those in 
Pearson et al. (2008). NMS was run in PC-ORD's "autopilot 
mode," which determined the optimal number of axes through 
a Monte Carlo test of 100 runs (50 runs on actual data, 50 
runs on random data). The final ordination included 99 runs, 
which were completed for the recommended number of axes. 
Environmental data were overlaid on the NMS plot to show 
correlation between lipids and geochemistry. Samples with no 
bGDGTs were removed from NMS and two-way cluster analy- 
ses to avoid the bias that may occur when a GDGT is assumed 
to be absent when, in fact, it may be present below detection 
concentrations. 

Non-parametric Mann- Whitney tests. Spearman's rho corre- 
lation coefficients, and select linear regressions were executed in 
IBM SPSS Statistics 19 to identify relationships among the GDGTs 
and geochemical analytes. Analyses were completed at the 0.05 
level of significance. 

RESULTS 

BRANCHED GDGTs IN HOT SPRING SEDIMENTS AND MATS AND 
ADJACENT SOILS 

bGDGTs were detected in 38 of 40 hot springs sampled, span- 
ning a temperature range of 31-95°C and a pH range of 6.8-10.7, 



across several distinct geothermal systems in the northwest Great 
Basin (Table 1; Figure 1; Table SI). Total bGDGT concentrations 
ranged up to 824ng/g dry mass in hot springs (Table 1); lower 
concentrations were observed in adjacent desert soils (Table SI; 
Figure 2). A strong, positive, linear relationship between log- 
transformed core and polar bGDGTs was observed (r^ = 0.692, 
sig < 0.001). Similar log-linear relationships existed for indi- 
vidual bGDGTs in hot spring samples, particularly bGDGT I, 
lb, Ic, II, and lib (Figure S2); however, no such relationship 
was clear for bGDGT III, while polar bGDGT lie, Illb, and 
IIIc were not detected in any hot spring samples. Samples from 
desert soils adjacent to the springs showed a similar relationship 
between log-transformed core and polar bGDGTs (r^ = 0.475, 
sig = 0.040); however, significant linear relationships for indi- 
vidual lipids were only obtained for bGDGT I and bGDGT lib 
(Figure S2). 

bGDGT I, lb, and Ic were present in many hot spring sam- 
ples at concentrations of one order of magnitude or more higher 
than in adjacent soils (Figure S2). A two-way cluster analysis 
based on the relative abundance of individual polar bGDGTs 
showed that hot spring samples had distinct bGDGT compo- 
sition compared with soil samples, with bGDGT I, lb, and 
Ic, in decreasing order, being the most common and abun- 
dant bGDGTs in hot spring samples and bGDGT II and III 
being the most common and abundant bGDGTs in soil samples 
(Figure 3). 

RELATIONSHIP BETWEEN bGDGT ABUNDANCE AND TEMPERATURE IN 
HOT SPRING SAMPLES 

bGDGTs were the dominant GDGTs in some samples, compris- 
ing up to 99.5% of total GDGTs (bGDGTs plus iCDGTs), whereas 
iGDGTs dominated other samples, with bGDGTs comprising 
< 1% (Figure 4). The relative abundance of bGDGTs and iGDGTs 
was related to temperature, with bGDGTs abundant at lower tem- 
peratures and iGDGTs dominating at higher temperatures (r^ = 
0.591, sig < 0.001; Figure 4). 

The absolute abundance of bGDGTs was also related to tem- 
perature, with most samples >70°C containing relatively low 
concentrations of core bGDGTs (<62 ng/g; Figure 5). This pat- 
tern was even clearer for polar bGDGTs, which are unstable in 
the environment, and therefore, generally interpreted to repre- 
sent living biomass (White et al, 1979). Polar bGDGTs were 
absent in all but one sample that was >70°C at the time of 
collection (Figures). A Mann- Whitney test confirmed the sig- 
nificance of 70° C as a cutoff for both core and polar bGDGT 
distribution within this dataset (core sig = 0.046; polar sig < 
0.001). Spearman's rank correlation analysis relating the abso- 
lute abundance of polar bGDGTs and geochemical measure- 
ments revealed significant negative relationships between polar 
bGDGT abundance and temperature (p = —0.582; sig < 0.001) 
and sulfide concentration (p = —0.405; sig = 0.014), and sig- 
nificant positive relationships between bGDGT abundance and 
dissolved oxygen (p = 0.508; sig = 0.001) and nitrate (p = 0.650; 
sig < 0.001) concentrations. Spearman's rank correlation anal- 
ysis with the absolute abundance of polar bGDGT I, lb, and 
Ic, analyzed separately, revealed similar relationships (Table S2). 
In contrast, the absolute abundance of other polar bGDGTs 
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FIGURE 2 I Absolute abundance of polar bGDGTs vs. core bGDGTs for 
soil and hot spring samples. Samples below the method detection limit 
for polar bGDGTs are not shown and were not used for in regression 
analyses (14 hot spring, 4 soil). Samples with <1 ng/g of polar or core 
iGDGTs, with negative loglO values are not shown; however, they were 
used for regression analyses (3 hot spring, 2 soil). 



had few significant correlations witli geocfiemical measurements 
(Table S2). 

The sample from the highest temperature in the Great Boiling 
Spring geothermal field, GBS-19, contained a low concentration 
of polar bGDGTs (2.2 ng/g dry mass), but the highest concen- 
trations of bGDGT lib (0.80 ng/g dry mass) and bGDGT III 
(0.80 ng/g dry mass) among hot spring samples. 

CO-VARIATION BETWEEN bGDGTs IN HOT SPRING SAMPLES 

bGDGT I, bGDGT lb, and bGDGT II, in order of decreasing 
abundance, were the dominant bGDGTs in hot spring samples 
(Table 1), comprising >70% of all polar bGDGTs in all sam- 
ples except for GBS-19. A two-way cluster analysis based on the 
relative abundance of individual polar bGDGTs was created to 
group geothermal springs and to examine the co-occurrence of 
individual bGDGTs in the environment (Figure 6). The cluster 
analysis along the horizontal axis revealed two major groups of 
samples. Group 1 included most low-temperature spring sam- 
ples, with all except one at <51°C. Group 2, in contrast, included 
sites ranging up to 70° C, with all except one collected from 
>50°C. Several other sites, mostly >70°C, had more unique 
bGDGT profiles and clustered separately. No geographical clus- 
tering was apparent, consistent with non-significant results in 
the Spearman's rank correlation analysis for the majority of 
chemical analytes, particularly redox-inactive species. The cluster 
analysis along the y-axis showed that the minimally methylated 
tetraethers, bGDGT I, bGDGT lb, and bGDGT Ic, co-varied, 
leaving a second, less tight cluster for the more highly methy- 
lated tetraethers bGDGT II, bGDGT lib, and bGDGT III. A 
heatmap associated with the two-way cluster analysis gener- 
ally supported both cluster analyses and allows visualization 
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FIGURE 3 I Two-way cluster analysis organizes hot spring samples 
(black) and soil samples (blue) according to similarity in polar bGDGT 
composition (horizontal dendrogram). bGDGTs are organized according 
to their presence and abundance in the samples (vertical dendrogram). 
Heat map colors indicate minimum (white) to maximum (red) abundance of 
bGDGTs, wherein colors are scaled by each lipid's maximum relative 
abundance in hot spring and soil samples. 



examined using nonparametric multidimensional scaling (NMS; 
Figure 7). In the NMS plot shown here, distance is proportional 
to multivariate dissimilarity in lipid composition (as calculated 
by a Sorensen distance matrix). Samples that are highly related 
are plotted closely together, while bGDGT types are plotted 
according to where the centroid of each lipid variable would 
occur. Physicochemical analytes are plotted as vectors, wherein 
the length and direction of a vector is proportional to its relative 
correlation with NMS axes. 

The 3-dimensional NMS model was a low-stress and statis- 
tically significant ordination (min stress = 5.921, p = 0.0196), 
which suggested strong underlying structure in lipid composi- 
tion. The NMS plot supported Group 1 and Group 2 identified 
by two-way cluster analysis and indicated a strong correla- 
tion between the lipid-derived NMS model and temperature 
(r^ = 0.438, correlation with NMS axis 3). Other strong corre- 
lations with the lipid-derived NMS model were oxygen (r^ = 
0.386, axis 3) and pH (r^ = 0.384, axis 3), both of which co- 
vary negatively with temperature in alkaline geothermal systems 
(Nordstrom et al., 2005). Individual bGDGTs were spread along 
a line originating near the center of the plot, populated by the 
minimally methylated bGDGT I, bGDGT lb, and bGDGT Ic, 
and terminating in the lower right quadrant with the highly 
methylated bGDGT III. The trend along axis 3 was in the same 
direction as the temperature vector and opposite to the pH and 
oxygen vectors, demonstrating a possible positive relationship 
between methylation and temperature and covariates of temper- 
ature; however, this trend was opposite the pattern of bGDGT I, 
bGDGT lb, and bGDGT Ic dominating in most hot spring sam- 
ples and may have been driven by GBS 19, which was an unusual 
sample that was dominated by bGDGT lib and bGDGT III. 
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FIGURE 4 I Percent bGDGTs (core plus polar) vs. temperature. Percent 
bGDGTs was calculated by dividing total bGDGTs by the sum of all GDGTs 
(iGDGTs + bGDGTs). 



of the distribution of each individual polar lipid within the 
dataset. 

Multivariate relationships between the relative abundance of 
individual, polar bGDGTs and physicochemical factors were 



DISCUSSION 

The abundance of bGDGTs in geothermal samples reported here 
is similar to those in geothermally heated soils at hot spring mar- 
gins ranging from 18 to 41°C (9-2,400 ng/g dry mass; Peterse 
et al, 2009) and in hot spring sediments in Great Boiling Spring 
(GBS), Nevada, ranging from 62 to 82°C (15-335 ng/g dry mass; 
Zhang et al., 2013) and Gulu and Yangbajing hot springs in Tibet, 
ranging from 52 to 83.6°C (up to 235 ng/g dry mass; He et al., 
2012). These bGDGT concentrations are also similar to those 
in non-thermal soils, such as those throughout the Pearl River 
watershed (10-770 ng/g soil; Zhang et al., 2012) but much lower 
than some terrestrial environments such as decaying peat (up to 
80 |xg/g dry mass; Weijers et al., 2006b). 

Our data here support the in situ production of bGDGTs 
in some Great Basin hot springs based on four lines of evi- 
dence: (1) the widespread distribution of bGDGTs in Great Basin 
geothermal springs, (2) the higher concentration of bGDGTs in 
many <70°C hot spring samples as compared with nearby desert 
soils, (3) the strong relationship between core and polar bGDGTs 
in geothermal samples (Liu et al., 2011), and (4) the distinctive- 
ness of soil and hot spring bGDGT profiles (Figure 3). These 
trends were particularly clear for bGDGT I, lb, and Ic, provid- 
ing strong evidence for their autochthonous production in the 
springs. However, some hot springs, particularly those >70°C 
had substantially lower abundance of bGDGTs than surrounding 
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FIGURE 5 I Absolute abundance of core and polar bGDGTs as a function 
of temperature and pH. Bubble areas (blue) are proportional to the absolute 
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represents the maximum absolute abundance in ng lipid g-1 dry mass. Gray 



circles represent zero abundance. Pink boxes delimit 70°C, at and below 
which bGDGTs are nearly ubiquitous in Great Basin hot springs. Statistical 
data refer to Mann-Whitney tests demonstrating the significantly higher 
abundance of bGDGTs <70°C. 



Relative Abundance (%) 

(normalized to maximum ot each lipid) 

min C 



Information Remaining (%) 
25 50 75 



100 



QQDQQQ 

X3 



Group 1 



1^ 



Group 2 



r BR(45) 
n- SV(60) 
DH(43) 
EVmat(44) 
SVX(50) 
SV(50) 
FGmat(40) 
SV(40) 
SSWof(51) 
DH(50) 
RHC2(62) 
SVX(60) 
DH{60) 
SSW/(40) 
FGlTBt(50) 
GBS(61) 
RHC1(52) 
V 

70) 

— GBSl9(95) 

I RHC3(69) 
^ SV(70) 



□□□□□□ 

□□□□□r 



□□□□□ 
□□□□C 

□□□□□L 

□□□□□r 

— UUL 



r RHC1(! 

__\ ' SSW(6 

|_r ssw(7 

l- S\/Xf7i 



■ DH(70) 
-■ ed( 



□□□□□ 

EVsed(44) 



Banc 

EDDDC 

nnnnc 

□BDC 



□□■■□L 
□□■□DD 

□□■□□r 

□□□■I 

□□■□c 



FIGURE 6 I Two-way cluster analysis organizes hot spring samples 
according to similarity in polar bGDGT composition (hiorizontal 
dendrogram). bGDGTs are organized according to their presence and 
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minimum (white) to maximum (red) abundance of bGDGTs, wherein colors 
are scaled by each lipid's maximum relative abundance in hot spring 
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soils, particularly polar bGDGTs (Figures 2, 5; Table SI). These 
observations are consistent with those reported in Schouten et al. 
(2007) and suggest that soil runoff may be a dominant source 
of bGDGTs in hot springs at temperatures >70°C. However, the 
predominance of autochthonous bGDGTs in lower temperature 
springs is consistent with evidence that bGDGTs can also be pro- 
duced in other aquatic systems such as lakes (Sinninghe Damste 
et al, 2009; Tierney et al, 2010; Sun et al., 2011; Wang et al. 
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FIGURE 7 I The NIVIS plot shows relationships among polar bGDGTs 
(relative abundance) and environmental variables. Samples (solid 
shapes) are arranged accordingly to similarity in polar bGDGT composition. 
Distance on the plot is proportional to dissimilarity in lipid composition. 
Lipids (red +) appear where the centroid of that variable would be placed 
within the NMS model. Environmental variables (gray) appear as vectors 
that indicate relative correlation with NMS axes. Colors of samples sites 
correspond to groups identified in the two-way cluster analysis. 



2012), rivers and estuaries (Zhu et al., 2011; Zhang et al., 2012), 
and geothermal springs (Zhang et al, 2013). 

All data point to temperature as a possible factor controlling 
bGDGT production and suggest that most bGDGT-producing 
thermophiles are moderate thermophiles rather than hyper- 
thermophiles. However, the effect of temperature on bGDGT 
production may be indirect, since incubation of lignocellulosic 
materials at 77 and 85°C in GBS led to in situ production 
of bGDGTs (Zhang et al, 2013). It is possible that bGDGT- 
producing thermophiles are members of consortia involved in the 
decomposition of complex biomass, which may only be found 
at high temperature on rare occasions, such as after delivery of 
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allochthonous plant material into the spring or during degrada- 
tion of photosynthetic mats after temperature up-shifts. The latter 
scenario could explain the variable presence of polar bGDGTs in 
>70°C sediments in GBS when it was sampled on different dates 
(Table 1; Zhang et al., 2013). Careful temperature monitoring and 
temporal sampling would be required to examine this possibility 
more rigorously. The temporal variability of the GBS geothermal 
field has been documented (Anderson, 1978) and is consistent 
with our own, more recent observations (e.g.. Cole et al, 2013; 
Peacock et al., 2013). Although polar bGDGTs were present at low 
concentrations in GBS 19, which has a stable source at 95°C, we 
urge caution interpreting that bGDGTs can be synthesized at that 
temperature since that site has a steep geothermal gradient and 
therefore it is difficult to infer the precise temperature at which 
lipids may have been produced. Furthermore, this is currently the 
only report of polar bGDGTs at temperatures >85°C and the pre- 
dominant lipids in that sample are heavily methylated, which is 
unusual in our dataset. 

A few other geochemical measurements correlated with abso- 
lute or relative bGDGT abundance: nitrate, sulfide, oxygen, and 
pH, in addition to temperature. Although a direct effect of these 
variables on bGDGT-producing bacteria cannot be ruled out, we 
caution against using these data to infer details about the spe- 
cific habitat and physiology of bGDGT-producing bacteria. These 
analytes co-vary with temperature in many geothermal outflows 
due to oxidation of reduced compounds present in the source 
water and gas (e.g., CO2) equilibration with the atmosphere (e.g., 
Nordstrom et al., 2005; Swingley et al., 2012) and their rela- 
tionship to bGDGT composition and abundance is difficult to 
infer. 

The absence of correlations between bGDGT abundance and 
any other geochemical analytes and the absence of biogeographic 
clusters in two-way cluster analysis and NMS plots (Figures 6, 7; 
Table S3) suggests that geological setting, and the resulting aque- 
ous and solid-phase geochemistry, are not primary drivers of 
bGDGT production within Great Basin hot springs. For exam- 
ple, no apparent clusters were formed by geographically distinct 
springs that are classified as Na-Cl springs (Great Boiling Spring, 
Sandy's Spring West, Rick's Hot Creek) nor Na-HCOj" -CI springs 
(Surprise Valley, Fly Geyser) (Anderson, 1978; Miller-Coleman 
et al., 2012). Our data are not consistent with a role of bGDGTs 
in stabilizing membranes in acidic conditions in these springs. In 
fact, bGDGT concentrations were directly proportional to pH, 
which is the opposite of what would be predicted if bGDGTs 
were serving a role to protect against low pH. This result does 
not rule out possible controls of geochemistry in bGDGT produc- 
tion in volcanicaUy active geothermal systems such as Yellowstone 
National Park or the Rehai Geothermal Field in Tengchong, 
China, each of which host springs with a much wider pH range 
than the Great Basin (Brock, 1978; Zhang et al, 1987). 

The wide range in the relative abundance of bGDGTs and 
iGDGTs in these springs (Figure 4) demonstrates great differ- 
ences in relative habitability for bGDGT-producing bacteria and 
iGDGT-producing archaea. The inverse relationship between 
bGDGT/iGDGT ratios and temperature (Figure 4) is in agree- 
ment with general trends favoring archaea at environmental 
extremes; an increase in the relative abundance of archaea 



along temperature transects within GBS has been documented 
previously (Cole et al., 2013). The dominance of minimally 
methylated bGDGT I, bGDGT lb, and bGDGT Ic among 
bGDGTs in this dataset is similar to the dominance of these 
three bGDGTs at higher temperatures in GBS (Zhang et al., 
2013) and Yellowstone hot springs (Schouten et al., 2007). Our 
results are also consistent with the high abundance of bGDGT 
II and III in geothermal soils in Surprise Valley and suggest 
that those bGDGTs may be produced by soil microorganisms 
(Peterse et al, 2009). Tibet hot springs had bGDGT III as 
the dominant bGDGT (He et al, 2012). These results sug- 
gest that factors other than temperature may have some influ- 
ence on the relative abundance of individual bGDGTs with 
different degree of methylation in geothermal ecosystems. The 
spreading of individual bGDGTs along axis 2 of the NMS 
plot (Figure 7), perpendicular to the temperature vector, also 
suggests factors other than temperature controlling degree of 
methylation. 

CONCLUSIONS 

The results presented here, along with previous studies (Schouten 
et al, 2007; Peterse et al, 2009; He et al., 2012; Zhang et al, 
2013) provide evidence for the production of bGDGTs by ther- 
mophilic bacteria inhabiting geothermal ecosystems, particularly 
bGDGT I, lb, and Ic. Although bGDGTs can be produced up to 
85°C (Zhang et al., 2013), the optimal conditions for bGDGT- 
producing microorganisms appears to be habitats <70°C, includ- 
ing both heated soils at hot spring margins and in sediments 
and mats within the springs themselves. As such, many bGDGT- 
producing bacteria appear to be moderate thermophiles, rather 
than hyperthermophiles. Temperature may be a master controller 
of the abundance of GDGTs in thermal springs, with moderate 
temperatures favoring a high relative abundance of bGDGTs and 
high temperature favoring high relative abundance of iGDGTs. 
However, the effect of temperature on bGDGT production may 
be indirect reflecting the presence of complex organic materials 
as incubation of lignocellulosic material at high temperatures led 
to production of bGDGTs at 77 and 85°C (Zhang et al, 2013). In 
contrast, geochemistry did not have a strong effect on bGDGT 
production, at least within the context of this study; however, 
other datasets, such as those from Tibetan hot springs, suggest 
global patterns of bGDGTs in geothermal systems may be more 
complex (He et al., 2012; Wang et al., unpublished data; Li et al., 
unpublished data). 
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